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ABSTRACT 

The condensations composing coronal rain, falling down along loop-like structures observed in cool 
chromospheric lines such as Ha and Ca II H, have long been a spectacular phenomenon of the solar 
corona. However, considered a peculiar sporadic phenomenon, it has not received much attention. 
This picture is rapidly changing due to recent high resolution observations with instruments such as 
Hinode /SOT, CRISP of SST and SDO. Furthermore, numerical simulations have shown that coronal 
rain is a loss of thermal equilibrium of loops linked to footpoint heating. This result has highlighted 
the importance that coronal rain can play in the field of coronal heating. In this work, we further stress 
the importance of coronal rain by showing the role it can play in the understanding of the coronal 
magnetic field topology. We analyze Hinode/ SOT observations in the Ca II H line of a loop in which 
coronal rain puts in evidence in-phase transverse oscillations of multiple strand-like structures. The 
periods, amplitudes, transverse velocities and phase velocities are calculated, allowing an estimation 
of the energy flux of the wave and the coronal magnetic field inside the loop through means of coronal 
seismology. We discuss the possible interpretations of the wave as either standing or propagating 
torsional Alfven or fast kink waves. An estimate of the plasma beta parameter of the condensations 
indicates a condition that may allow the often observed separation and elongation processes of the 
condensations. We also show that the wave pressure from the transverse wave can be responsible for 
the observed low downward acceleration of coronal rain. 

Subject headings: magnetohydrodynamics (MHD) - Sun: corona - Sun: flares - waves 



1. INTRODUCTION 

There is now increasing evidence that active regions 
in the Sun have their heating concentrated mostly at 
lower atmospheric regions, from the lower chromosphere 
to the lower corona. The excess densities found in most 
observed coronal structures such as coronal loops put 
them out of hydrostatic equilibriur n, a state that ca n 
be explained by f ootpoint heating ('Asch wandenI l2QQl[ ). 
iHara et all (|2QQ8[ ) using the Hinode /EIS instrument, 
have shown that active region loops exhibit upflow mo- 
tions and enhanced nonthermal velocities at their foot- 
points. The characteristi c over-density has also b e en de- 
duced seismologically by I Van Doorsselaere et al.l (|2QQ7[ ) 
by studying the fundamental to the second harmonic 
period ratio P1/P2 of standing transverse oscillations 
in loops. Recently, iDe Pontieu et al.l (|2Qllf ) have high- 
lighted the importance of the Hnk between the photo- 
sphere and the corona by showing that a considerable 
part of the hot coronal plasma could be heated at low 
spicular heights, thus explaining the fading character 
of the ubiquitous 'type II spicules'. Further evidence 
of footpoint heating is put forward by the presence of 
cool structures in the active region coronae, such as fila- 
ments/prominences or coronal rain, two phenomena that 
may share the same formation mechanism but which 
seem to differ on the structure of their underlying mag- 
netic field topology, leading to different observational as- 
pects such as dynamics, shapes and hfetimes. 

Both, prominences (or filaments if observed on disk 
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rather than at the limb) and coronal rain correspond 
to cool and dense plasma observed at coronal heights 
in chromospheric lines such as Ha and Ca II H and K. 
But while the plasma in prominences is suspended in the 
corona against gravity making the structures long-lived 
(days to weeks), coronal rain is observed falling down 
in tim escales of minutes (|Schrijverll2QQll : iDe Groof et al.l 
I2QQ4I ) along curved loop-like trajectories. The mechan- 
ical stability and thermodynamic properties of promi- 
nences are linked with the underlying magnetic field 
topology, and thus the main difference between coronal 
rain and prominences could be a difference in coronal 
magnetic field configuration. This question awaits fur- 
ther proper investigation. 

While prominences have been studied extensively in 
solar physics, few observational studies exis t of coronal 
rain s ince their di scovery in the early 1970s (jKawaguchil 
I197QI : lLerovlll972f ). This lead to the belief that coro- 
nal rain is a rather uncommon phenomenon in active 
region coronae. Furthermore, coronal rain is often erro- 
neously attributed to prominence material falling back 
from coronal heights following a prominence eruption. 
However, recent high resolution observations with in- 
struments such as CRISP of SST, SOT of Hinode or 
SDO reveal coronal rain to be dynamic, short-lived (1-10 
minutes), small sized (200 km or less) strand- like struc- 
tures that are ubiquitously pres ent over active regions 
(jAntolin et al.l 12011. in preparati on) . The detection of 
coronal rain therefore requires high spatial and temporal 
resolution observations that have only recently become 
available. The same type of observations have shown that 
prominences are composed of a myriad of fine threads, 
outlining a fine-scale structure of the magnetic field, and 
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the presence of flows along the threads (iHeinzel fc Anzerl 
^QQllLin et aI]l2(M[MM[Ur][m^ et aL 2008). 

The frequency of cor onal rain has profoun d implications 
for coronal heating (j Antolin et al.l l2010f ) . In order to 
quantify the true occurrence frequency and importance 
of coronal rain in active region loops we will require fur- 
ther effort gathering a statistically signiflcant number of 
high resolution observations in different wavelengths. 

Numerical simulations have shown that coronal rain 
and prominences are most likely the result of a phe- 
nomenon of thermal in sta bility , also known as 'catas- 
trophic cool i ng' (Hildner 19 741: lAntiochos et al 
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Antolin et al.l l2010f ). Loops with footpoint heating 
present high coronal densities and thermal conduction 
turns out to be insufficient to maintain a steady heat- 
ing per unit mass, leading to a gradual decrease of the 
coronal temperature. Eventually recombination of atoms 
takes place and temperature decreases to chromospheric 
values abruptly in a timescale of minutes locally in the 
corona. This is accompanied by local pressure losses 
leading to the formation of condensations which become 
bright or dark if observed towards the limb or on disk, 
r espectively. 

lAntolin et al.l ()2010f ) showed that Alfven wave heat- 
ing, a strong coronal heating candidate, is not a pre- 
dominant heating mechanism in loops with coronal rain. 
When propagating from the photosphere into the corona, 
Alfven waves can nonlinearly convert to longitudinal 
modes through mode conversion due to density ffuctu- 
ations, wave-to- wave interaction, and deformation of the 
wave shape during propagation. These modes subse- 
quently steep i nto shocks and heat the plasma uniformly 
along the loop (Morivasu et al. 2004; " Antolin fc Shibatal 
1201 Ol : IVas heghani Farahani et al. 20 ll]), thus avoiding 
the loss of thermal equilibrium in the corona. Coro- 
nal rain is often observed falling down at speeds much 
lower than free fall speeds resulting from the effective 
gravity along loops ('Schriiver' 200TI:|pe Groof et al.ll2004l: 
lie Groof et al. 2005; Miiller et al.l 120051 : lAntolin et all 
[2Q10i) . Simulations have shown that the effects of gas 
and ma gnetic pressure may explain the observed dy- 
namics (iMackav fc Galsgaar dl [200TI : iMiiller et al.l [20031 : 
lAntolin et al.l I2010D . Here we show that the observed 
wave pressure from a transverse wave may also account 
for decreased accelerations. 

Magnetohydrodynamic (MH D) waves are fre q uently 
observed in prominences ((Ramsey fc SmithI 119661: 




Oliver Ballesterl [200l iFoullon etal.l 120041 
120091 : lOkamoto etHI 1200?!; If err 
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and coronal loop struc t ures (lA schwanden et al. 
iDe Moortel et al.l [200(1 IVan D oorsselaere et al. 
lErdelvi TarovanI 120081 : IVerwichte et all 
leading to the determination of the in- 
ternal physical conditions through the develop- 
ment of analy tical th eory and numerical m o deling 
([Roberts et al.f 119841: iNakariakov Verwichtd 120051: 



[Ballester '2006; Andries e t al.l l2009riTarovan fc Erdelvil 
[2009; Arregui et al. 2010), a technique dubbed coronal 
seismology. The determination of the physical prop- 
erties of the corona through which MHD waves travel 
depends on the correct interpretation of the observed 
signatures, correct identiflcation of the wave mode and 
an MHD wave model that provides robust seismological 



measurements. For waves in transients objects such as 
spicules and fllament flbres, the role of a wave guide has 
b een debated. 

lOkamoto et al.l (|2007[ ) analyzed transverse oscillations 
running through prominence threads observed by Hin- 
ode /SOT in the Ca U H-line at the limb of the Sun. The 
reported mean periods for the waves are between 130 s 
and 240 s, (horizontal) oscillation amplitudes between 
400 km and 1770 km, transverse (vertical) velocities be- 
tween 5 km s~^ and 15 km s~^ and an estimated wave 
speed larger than 1050 km s~^ leading to a magnetic fleld 
of 50 G in the prominence. Mini mum Alfven speeds in 
the prominence were estimated by iTerradas et al.l ()2008l ) 
to be between 120 km s~^ and 350 km s~^, depend- 
ing on the local magnetic fleld, the total lengths of 
the magnetic fleld lines in the prominence, and the ra- 
tio between th e loca l and external (coronal) density. 
lOkamoto et al.l ()2007[ ) flrst interpreted the oscillations 
as Al fven waves running th roug h the prominence. How- 
ever, [Tbrradas et al.l (|2008[ ) and [Van Doorsselaere et al.l 
(|2008aj ) have argued that the only solution among fast 
waves that gives rise to a displacement of a magnetic flux 
tube a xis is the kink mode. Furthermore. ITerradas et all 
(|2008[ ) have shown that the periods of the kink mode are 
rather insensitive to the presence of steady flows along 
t he threads. 

lOfman fc Wand (|2008l ) have analyzed an event with 
Hinode/ SOT in the Ca ll H line, similar to the one in 
the present work. Transverse oscillations are observed 
in a loop with flows. In this case, the cool material is 
ejected at very high speeds (74 — 123 km s~^) from one 
footpoint to the other, and is related to a flare happen- 
ing close-by, which may also be the cause for the oscilla- 
tions. The waves are interpreted mostly as fundamental 
modes of standing kink oscillations, although some of 
the observed threads display dynamics more consistent 
with propagating fast magnetoacoustic waves. Coronal 
seismology is performed assuming a density in the range 
(1 — 5) X 10^ cm~^, leading to coronal magnetic flelds of 
20 ± 7 G. The analyzed loop does not seem to be subject 
to catastrophic cooling during the observed time, and 
thus the flow is of a different nature than that of the 
present work. The cause for the observed oscillations in 
our case seems to be different as well, since no energetic 
phenomenon is observed. 

In this wo rk we analyze the sam e high resolution ob- 
servations of lOkamoto et al.l ()2007l ). but concentrate on 
active region loops in the foreground, unco nnected to the 
promi nence, that exhibit coronal rain. In lAntolin et al.l 
(|2010[ ) the observational analysis concentrated on loops 
to the north of the visible sunspot. Here, we will focus 
on one loop on the south side of this sunspot and which 
exhibits a peculiar phenomenon. We present the flrst ob- 
servational analysis of transverse oscillations of threads 
in a loop subject to coronal rain. The paper is orga- 
nized as follows. In Section [2] we describe the data set of 
Hinode/ SOT ^ present statistics of velocities and accelera- 
tions for the falling coronal rain and analyze the observed 
oscillations in the loop. In Section [3] we proceed to dis- 
cuss the observational results, giving interpretations for 
the wave properties and their nature, and flnalize in Sec- 
tion [4] with the conclusions of the work. 
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Figure 1. Active Region NOAA 10921 on the west limb observed by Hinode/SOT in the Ca II H band on the 9 November 2006 between 
19:33 and 20:44 UT. The curves denote some of the paths traced by coronal rain. The soh d curves conform the loop studied in the present 
work, while the loops outlined by the dashed curves were studied in lAntolin et al.l (|2010I V The dotted curves mark the presence of other 
loops which may be interacting with the loop studied here. 
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Figure 2. Active Region NOAA 10921 observed by Hinode/XRT with the Al Poly filter on the 9 November 2006 at 19:59 UT, roughly 
half an hour before the observed coronal rain in the studied coronal loop. The solid and dashed lines mark the position of the loop, as in 
Figure [3] 
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Figure 3. Loop with coronal rain observed by Hinode/SOT on 
the 9 November 2006 above AR 10921. 



2. OBSERVATIONS OF CORONAL RAIN WITH 
HINODE/SOT 

2.1. Velocities and accelerations 

The observatio ns with the Solar Optical Telescope 
(SOT) of Hinode (jTsuneta et al.ll2QQ8[ ) are in the Ca II H 
band, on 9 November 2006 from 19:33 to 20:44 UT with a 
cadence of 15 s and a spatial resolution of 1.22X/D 0^2, 
and focused on NOAA AR 10921 on the west limb. A 
variance of the images over a part of the time interval 
is shown in Figure [TJ This data set has become fa- 
mous among solar limb observations. The set shows the 
presence of an active region promine nce which e xhibits 
interesting oscillatory behaviour fOk amoto et al.| [2Q07). 
Figure [2] shows a Hinode /XKF observation of the same 
region at 19:59 UT with the Al Poly filter. The square 
in the Figure corresponds to the SOT field of view. 

Additionally, on the foreground of the prominence var- 
ious loops exhibiting coronal rain have been observed. A 
statistical study of the loops outlined in dashed curves 
in Figure [1] can be found in Antolin et al. (2010). The 
observed loops are located north of the sunspot, have 
lengths between 60 Mm and 100 Mm, and exhibit coronal 
rain continuously. The condensations composing coronal 
rain display a broad distribution of velocities (between 20 
km s~^ and 120 km s~^) that put in evidence both accel- 



eration and deceleration processes in the loops. The av- 
erage accelerations were found to be lower than that pro- 
duced by gravity, indicating the presence of other forces, 
possibly of magnetic origin. 

The focus of our paper is set on a coronal loop lo- 
cated south of the sunspot, outlined in solid curves in 
Figure [TJ which can be observed in the Ca II H band of 
Hinode/SOT thanks to the coronal rain occurring in the 
loop. Figure [3] shows the subset of the entire field of view 
corresponding to this loop. The loop is visible for about 
half an hour towards the end of the observation set. We 
have plotted the variance of the image over the period 
of time it becomes visible. Assuming that the geometry 
of the loop is close to that of a semi-torus, we see from 
Figure [3] that the plane of the loop makes a significant 
angle with the plane of the sky, being roughly directed 
along the line of sight and that it is slightly inclined with 
respect to the vertical. The coronal rain can be observed 
basically from the apex of the loop, located 25 ± 5 Mm 
above the surface, leading to a loop length of 80 ± 15 Mm 
assuming a circular axis for the loop. In the Hinode /XBT 
image of Figure [2] we have outlined the position of this 
loop. 

We have tracked down the condensations along the 
loop with the help of CRISPEX (CRisp SPectral EX- 
plorer) and Timeslice ANAlysis Tool (TANAT) 0, two 
widget based tools programmed in the Interactive Data 
Language (IDL), which enable the easy browsing of the 
image (and if present, also spectral) data, the determi- 
nation of loop paths, extraction and further analysis of 
length- time diagrams. 

The condensations that can be tracked from high up 
in the corona down to chromospheric heights normally 
have a considerable thickness of about half a megameter. 
However, separation and elongation of the condensations 
generally occurs during the fall, leading to very thin and 
elongated condensations tracing strand-like structures. 
Many smaller condensations can be observed at various 
heights but are normally too faint to be followed clearly 
along their paths. In the observed loop a total of 28 con- 
densations can be easily tracked to the chromosphere. 
Velocities and accelerations are derived with the help of 
length- time diagrams, which are shown in the left and 
middle panels of Figure HI Since the velocity of the 
condensations varies along their paths, when possible, 
multiple velocity measurements at different heights are 
made, which allows to estimate the acceleration at dif- 
ferent heights. The statistics in this loop are similar to 
those of the other loops on the north side of the sunspot. 
A broad distribution of velocities between 20 km s~^ and 
100 km s~^ with a mean below 40 km s~^ is obtained. 
The accelerations have in average lower values and are 
concentrated around a mean of 0.056 km s~^. Now, the 
change of the average effective gravity for a blob in free 
fall from the top of an ellipse with respect to its ellipticity 

can be calculated easily as {geff) = f Jq^'^ 9q cos 0{s)ds, 
where 0{s) is the angle between the vertical and the tan- 
gent to the path and s is a variable parametrizing the 
path. It is found that for a ratio of loop height to half 
basehne between 0.5 and 2, (geff) varies roughly be- 

^ The actual code and further information can be found at 
|http: //www, a stro, uio.no / ~gregal / crispex / index, html 
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Table 1 

Oscillation features for the events 
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Velocity [km/s] 
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Figure 4. Histograms of velocity (left), acceleration (middle) and 
(projected) height versus velocity (right) for the coronal rain ob- 
served with Hinode /SOT. 



tween 0.132 km s~^ and 0.21 km s~^, values that are 
significantly larger than the observed average value. A 
few cases of faster acceleration, as well as decelerations 
and constant falling velocities are also observed. In the 
right panel of Figure |4] we have plotted the heights of the 
measurements for the condensations with respect to their 
velocity at that height. The solid and dashed lines in the 
figure correspond to the free fall speed under the action 
of the solar surface gravity and with the mean observed 
acceleration, respectively. 

Since Doppler velocities are not available in the cases 
analyzed in this paper we can only measure projected 



Events Period 



Amplitude 
[km] 



Transversal 
velocity [km/s] 



112 ± 29 
171 ± 11 
118 ± 15 
143 ± 17 

165 ± - 
198 ± 49 
176 ± 23 

84 ± 8 



245 ± 148 
515 ± 135 
308 ± 136 
324 ± 125 
351 ± 55 
371 ± 205 
369 ± 255 
305 ± 119 



4.5 ± 2.5 

6.1 ± 1.4 

5.3 ± 1.3 

4.4 ± 1.2 

4.2 ± 0.1 

3.6 ± 2.0 
3.8 ± 2.7 
7.6 ± 4.0 



Note. — Periods, amplitudes and transversal ve- 
locities with respective standard deviations for the 
8 detected oscillatory events. Since event 5 is only 
observed to oscillate through 1 period before fading 
out, the calculation of the standard deviation is not 
possible for this case. 

velocities without further assumptions on the geometry 
of the loop. The velocities and accelerations in the pan- 
els are thus lower estimates of the true values. We can 
however make an estimate of the errors. From Figure [3] 
the projected distance on the plane of the sky between 
the two footpoints of the loop is estimated to be / 12 
Mm. This implies an angle between the line of sight and 
the plane of the loop of roughly 14°. Naming h and H 
the height of the measurement and the total height of 
the loop, respectively, the obtained error for each mea- 
surement results 



error : 



hl/2H 



V^l - {hl/2HY 



(1) 



Since the upper and lower heights of the measurements 
are around 17 Mm and 2 Mm respectively. Equation ([1]) 
gives an error of 16.5% and 2% for the upper and lower 
velocity measurements, respectively. This results in an 
error = 12.5% for the acceleration. Hence the true mean 
acceleration is roughly 0.049 km s~^. 

2.2. Oscillations 

As the condensations fall, separation and elongation 
processes occur, which results in several strand-like struc- 
tures resolved in the loop. The strands are observed 
to oscillate transversally. Figure [5] shows the time slice 
along the outlined loop of Figure [3l where the transverse 
length refers to the perpendicular distance to the dot- 
ted line in Figure [3l from dashed line to dashed line. 8 
oscillation patterns can be clearly observed, over which 
we plot in color the (projected) distance where they hap- 
pen from the apex of the loop. There is a clear in-phase 
oscillation for the strands 1, 2, 3, 4 and 6 between the 
time period [15, 22] min in the figure. These strands 
become all visible simultaneously at roughly 8 Mm from 
the apex. 

In Table [T] the estimated periods, peak to peak ampli- 
tudes and transversal displacement velocities of the oscil- 
lations are shown. The periods lie between 100 s and 200 
s, amplitudes are all roughly below 500 km and transver- 
sal velocities are between 4 km s~^ and 8 km s~^. The 
standard deviations are calculated by taking into account 
the widths of the strand-like condensations and the pos- 
sible errors that they involve. Since these widths can be 
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Figure 5. Time slice across the loop. The transversal width corresponds to the distance between the two dashed lines in Figure [S] along 
a perpendicular to the dotted line. The time interval includes the time when coronal rain is observed. 8 oscillations can be detected. We 
repeat the figure, plotting in color over the oscillations the length from the apex where they occur. 



up to 500 km, this leads to large uncertainties in the oscil- 
lation amplitude measurements. This is also reflected in 
the calculation of the distance from the apex where the 
oscillation occurs (plotted in color in Figure [5]), where 
the distance does not always vary smoothly. 

Along their paths from the apex towards the chromo- 
sphere, the strand-like condensations reach a maximum 
separation between each other at an apparent distance 
of 4 Mm to 8 Mm from the apex of the loop, after which 
they gradually converge in the lower part of the loop's 
leg to a common footpoint about 1 Mm wide in the chro- 
mosphere. The maximum observed separation between 
the strands is close to 5 Mm wide. If the condensations 
indeed follow the magnetic field, this implies a magnetic 
geometry with a cross-sectional area expansion factor of 
at least 25 in 20 Mm height between chromosphere and 



corona. 

In the present observations by Hinode/SOT several 
strand-like structures in a loop are observed to oscillate 
in-phase (synchronous). This points to either one trans- 
verse magnetohydrodynamic wave that affects all blobs 
as part of one monolithic loop, or either multiple trans- 
verse waves excited in separate strands but excited in 
phase by a common large-scale source. Possibly, the os- 
cillations are not only confined to the specific structures 
but can involve the larger coronal region, and thus could 
be related to the prominence osci llations visible in the 
background (jOkamoto et al.l [20071 ). However, although 
the oscillation periods are similar, the difference in mode 
polarisation and the absence of oscillations in loops at 
the north side of the sunspot imply that the line-of-sight 
distance to the prominence may be too large to expect a 
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Table 2 

Wave properties and magnetic field 



10 15 20 

Position [Mm] 



25 



30 



Figure 6. Oscillation amplitude (transversal displacement from 
peak to peak oscillation) with respect to position along the loop. 
The position is calculated from the projected length from the apex 
(0 is at the apex), assuming a circular loop of height 25 Mm. The 
calculated points correspond to averages over boxed regions around 
each peak to account for the width of the strand. The error bars in 
position correspond to the standard deviations. The error bars in 
amplitude correspond to the width of the strand at each oscillation 
peak. The solid line corresponds to a fitted to the data of the first 
harmonic sine profile. 



common excitor. 

Information about the longitudinal and propagatory 
nature of the wave producing the observed oscillations 
can be obtained by analyzing the change of the oscillation 
amplitudes with respect to the position along the loop. 
This is plotted in Figure [6l Each point corresponds to 
a peak-to-peak amplitude of a specific oscillating strand 
and the corresponding position along the loop where it 
happens, corrected for projection effects assuming a cir- 
cular loop of 25 Mm in height. The error bars in position 
correspond to the standard deviation of all the calculated 
positions in a boxed region around each peak. The er- 
ror bars in amplitude are set equal to the width of the 
strand at each oscillation peak, which explains why they 
are large. 

Since the oscillations in the loop can only be observed 
when the condensations are falling it is difficult to di- 
rectly ascertain whether the agent is a propagating or a 
standing wave. Figure [6] shows that the condensations do 
not appear to be oscillating when they are at the apex 
of the loop, nor in the lower part of the loop close to the 
footpoint. Furthermore, the amplitudes indicate a max- 
imum at roughly half way along the loop leg (one fourth 
of the total loop length), which correspond to signatures 
of the first harmonic of a standing mode. For compar- 
ison, the solid line in the figure corresponds to a fitted 
sine profile to the data, which is the profile that a first 
harmonic would have. Alternatively, one could envisage 
a propagating wave packet, propagating up or down, the 
maximum amplitude meeting the condensations half way 
through the loop's leg. However, it is difficult to see how 
a propagating wave may reach a maximum amplitude at 
a given height independent of wave amplitude. There- 
fore, this scenario is less likely. 

Apart from creating a transversal oscillation of the con- 
densations in the loop, the observed waves may also have 
an effect on the general dynamics of coronal rain. As 



Events 



Phase speed 
[km s-i] 



Magnetic field 

[G] (Pe = PO) 



Magnetic field 

[G] (Pe = 0) 



761 ± 264 
461 ± 30 
677 ± 89 
555 ± 67 
476 ± - 
411 ± 91 
452 ± 60 
942 ± 96 



21.5 ± 7.5 
13 ± 0.8 
19 ± 2.5 

15.7 ± 1.9 
13.5 ± - 

11.6 ± 2.5 

12.8 ± 1.7 

26.7 ± 2.7 



15.2 ± 5.3 
9.2 ± 0.6 
13.5 ± 1.8 
11.1 ± 1.3 

9.5 ± - 
8.2 ± 1.9 

9 ± 1.2 
18.8 ± 1.9 



Note. — Average phase speeds, energy fiuxes and mag- 
netic fields for each oscillation event in Table [U with the 
corresponding standard deviations. The magnetic fields 
are calculated according to Equation [2] for two limiting 
cases: an exterior to interior density ratio of 1 {third col- 
umn), or a density ratio of (fourth column). 



seen in the middle panel of Figure |4] the observed ac- 
celerations do not have a broad distribution but concen- 
trate around a low mean value of 0.056 km s~^. The 
range of the distribution is considerably smaller than 
that of the coronal rain observed in the coronal loops 
on th e north side of the SOT field of view (|Antolin et al.l 
l2QlQf ). This may be a signature of a specific force acting 
in the upward direction. A net upward wave pressure 
should be present, which could be both in the (upward) 
propagating wave scenario or in the standing wave sce- 
nario. In the later, a first harmonic would produce in 
the upper first half of the leg a downward acceleration, 
followed by a deceleration in the lower part of the leg, 
which is the portion of the loop that is mostly observed. 
A transverse magnetohydrodynamic wave would exert 
an average acceleration on the plasma proportional to 
(87rp)~^A5^/A/i, where B± is the transversal compo- 
nent of the magnetic field, p is the density of the plasma 
and h refers to a particular height. The average value 
for the effective gravity along the loop that a conden- 
sation would feel is ^Qq/tt = 0.174 km s~^, assuming 
a circular loop, which implies an observed average de- 
celeration of 0.118 km s~^. In order to obtain the 
observed deceleration purely by means of the wave pres- 
sure, assuming a typically high coronal number density 
of 3 X 10^ cm~^ for active region loo ps subject to catas- 
trophic cooling (|Antolin et al.|[2QTQl ), we would need a 
variation of AB± 0.4 G in a 1 Mm height difference. 
1.5 dimensional MHD simulations of both, standing tor- 
sional Alfven waves heating coronal loops, and propa- 
gating Alfven waves powering the solar wind, typically 
exhibit suc h gradients in the azimuthal magnetic field 
comp onent (jAntolin fc S hibata'l 2010l : ISuzuki fc Inutsukal 
120061 ). Furt hermore, an alytical analy sis and numerical 
modeling bv lTerradas fc OfmanI (|2004f ) have shown that 
the ponderomotive force from MHD waves (specifically 
from standing waves) can create fiows and lead to sig- 
nificant density enhancements at locations of maximum 
amplitude. 

3. DISCUSSION 

During the time interval in which the blobs are ob- 
served no phase shift is detected in the oscillations, re- 
inforcing the standing wave scenario. In case that prop- 
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Figure 7. Tracking of condensations 1 (upper) and 2 (lower) of Figure O The length in the x axis corresponds to the (projected) length 
from the apex of the loop. The length in the y axis corresponds to the same transversal length as in Figure O The times (in minutes) are 
set in the low left corner of each panel. Notice that each condensation does not oscillate up and down as a single structure but exhibits a 
different oscillation amplitude along its length. 



agating waves are playing a role in the oscillation of the 
blobs, these have to propagate with a phase speed larger 
than that set by the uncertainty in the measurements. 
Taking for instance the case of strand 1, which is ob- 
served over a distance of about 12 Mm according to Fig- 
ure O with an uncertainty of 1 time step we have a lower 
bound for the phase speed on the order of 800 km s~^. 
For the onset of a first harmonic in the loop, the to- 
tal wavelength is equal to the loop length, leading to 
phase speeds between 400 km s~^ and 1000 km s~^ (see 
Table [J) • To give a measure of the energy contents of 
these waves, we calculate the wave energy flux for the 
case of strand 1 as would be if the wave was propagating 
with half the amplitude. For this case, the phase speed 
is Vph = 760 ± 264 km s~^, giving an Alfven speed of 
va = Vph/'^^^'^ = 540 ± 186 km s~-^. The energy flux is 
then given by £^flux = ^P'^t '^a, where p = finemp is the 
average coronal density of the loop (taking an electron 
number density of rig = 3 x 10^ cm~^ and /i = 1.27 the 
effective particle mass with respect to the proton mass 
mp), and Vt is the transversal velocity calculated in Ta- 
ble [H Replacing with the numerical values we obtain an 
energy flux of (4.9 =b 5.9) x 10^ erg cm~^ s~^, where the 
large uncertainties are due to the large standard devia- 
tions in the observed periods and transversal velocities. 



The calculated value is below the estimated necessary 
energy of 10 ^ erg cm~^ s~^ for heat ing the active re- 
gion corona (|Withbroe fc Noveslll977[ ). In order to as- 
sess how much energy is transferred from the standing 
wave to the plasma, information about the wave damp- 
ing is necessary. However, in our case, it is difficult to 
find conclusive evidence for damping of the oscillation 
because in tracking a falling condensation the informa- 
tion of the height-dependence of the oscillation profile 
and the time-evolution are intertwined. 

Since we are observing the loop in a state of cooling, 
long after the heating has taken place, it is possible that 
the initial oscillations were faster, thus implying a much 
larger energy fiux for the initial propagating waves set- 
ting up the standing oscillations. Let us assume this 
scenario, and a sufficient energy flux to heat the corona 
of 10^ erg cm~^ s~^ at the beginning of the 'conden- 
sation' cycle (the cycle where the catastrophic cooling 
takes place, also known as 'limit cycle'). An initial coro- 
nal electron number density of 10^ cm~^ (a rough av- 
erage for active region loops) and a coronal magnetic 
fleld of 5-, 10 G (tak ing Equation [2l see also Figure 12 
in „Nakariakov fc Ve rwichte 2005) leads to a transversal 
velocity for the oscillations of at least 22 km s~^. This 
means that in a time interval of about 20 to 30 min- 
utes (which is an estimation of the 'condensing' phase 
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Figure 8. Sketch illustrating the possible interpretations to the oscillatory events. First, a standing kink mode corresponding to a 
first harmonic, that displaces the axis of a single monolithic wave guide occupying the whole magnetic structure (left). Second, multiple 
standing kink modes, each within individual strands in the magnetic structure guided by local density enhancements associated with each 
blob (middle). Third, a torsional Alfven wave that occupies one or several flux surfaces within the single monolithic loop producing a 
swaying of the individual strands from the twisting motions, without displacing the main axis of the loop (right). The blobs in red represent 
the condensations falling along the loop, from the apex towards the right footpoint. 



time prior to catas trophic cooling, based on results of 
lAntolin et all \201(t ) our oscillations have slowed down 
to about 1/3 to 1/4 their initial values at the begin- 
ning of the cycle. During the catastrophic cooling we 
do not observe any damping (mostly due to the fact 
that we cannot observe the same portion of the loop 
over a significant time interval), but it is possible that 
the main damping has occurred already during the con- 
densing phase. As decaying mechanisms, resonant ab- 
sorption has been shown to be an effective mech anism 
(|Ruderman fc RobertsI [2QQ2I : iGoossens et all 120021 ). As- 
suming an exponentially decaying amplitude, and taking 
an observed mean transversal velocity of 5 km s~^ (ac- 
cording to Table [Tj), we obtain a damping rate of (0.05 - 
0.075) min~^ for a time interval of (20 - 30) min. Thus 
leading to a decaying time of (13.5 - 20) min, which is 
in the order of whathas been reported (matching the pe- 
riods observed here) in simila r loops (Nakariakov et al. 
119991 : lAschwanden et al.l l2002^. This calculation is how- 
ever strongly dependent on the time of the condensation 
cycles (limit cycles), which are still subject to debate. To 
confirm such scenarios direct evidence from EUV imagers 
in future events will be required. On the other hand, 
if the initial energy of the wave is not dissipated but 
conserved, it is worth noting that the transversal veloc- 
ities of the oscillation are not altered significantly by an 
increase of density through the condensation processes, 
since it can be seen from conservation of wave energy fiux 
that the transversal velocity amplitude Vt depends only 
weakly on density increases as p as Vt ex 

The determination of the local magnetic field seis- 
mologically depends on the correct identification of the 
wave. This, in turn, depends on the existence of a waveg- 
uide in the loop. Prior to a catastrophic cooling event, 
simulations show that the loop undergoes a phase of in- 
creasing coro nal density and slowly decreasing co ronal 
temperature (lAntiochos et all 119991 : iMiiller et aH 120041 : 
lAntolin et al.ll2010f ). Figure [2| which corresponds to an 



XRT observation roughly half an hour before the ob- 
served coronal rain, shows that most of the region in the 
SOT field of view is composed by plasma with temper- 
atures above a million degrees. Our loop is not directly 
visible, which may be due to the coohng phase prior to 
the catastrophic cooling event. During this phase, which 
constitutes most of the cycle, the average coronal density 
reaches values that are considerably higher than that of 
the average exterior corona. According to simulations, 
this dense loop state, which can constitute a waveguide, 
is maintained for times on the order of tens of minutes, 
enough for the onset of standing modes from fast mag- 
netoacoustic waves. In our case, the observed transverse 
displacement can be associated with one of three types 
of waves, sketched in Figure [HI First, the wave is a kink 
mode that displaces the axis of a single monolithic wave 
guide that occupies the whole magnetic structure. Sec- 
ond, there may be multiple kink modes, each within indi- 
vidual strands in the magnetic structure guided by local 
density enhancements associated with each blob. Third, 
the wave is a torsional Alfven wave that occupies one or 
several fiux surfaces within the single monolithic loop. 
The main axis of the loop is then not displaced but we 
observe the swaying of the strands from the twists or 
torsional motions within fiux surfaces. A slow magne- 
toacoustic wave is excluded from consideration because 
it is essentially a longitudinal mode. 

For all cases the average coronal magnetic field in the 
loop can be determined from: 



Bo 



2-K- 



L 



(2) 



(formula 31 in iNakariakov fc Verwicht^ I2005L corrected 
for the first harmonic), where pe and po ^ire the exte- 
rior and interior densities respectively. This equation is 
valid under the assumption that the width of the loop 
is much smaller than the total length, which holds ac- 
cording to the estimates of section [2?2l The case of pe=^ 
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corresponds to the lower limiting case for a kink mode, 
whilst pe=po would correspond to the case of a torsional 
Alfven wave. In Table [2] we have estimated the magnetic 
field in both cases, assuming the height of the loop to be 
25 Mm and a loop number density of 3 x 10^ cm~^, a nor- 
mal value of dense loops subject to catastrophic cooling 
(o ur loop here is not s o different from the modeled loop 
in I Antolin et alJ |2Q1Q| ) . Taking the two limiting cases 
of exterior-to-inside density ratios of (kink mode sce- 
nario) and 1 (torsional Alfven wave scenario) we obtain 
an average coronal magnetic field of 11.5 G and 16.3 G, 
respectively. Importantly, when applying Equation [2] we 
make the assumption that the blobs themselves do not 
carry sufficient inertia to affect the oscillation itself. A 
scenario in which this assumption does not hold is worth 
investigating, and will be addressed in a subsequent pa- 
per. 

It is interesting to note that the collective behavior of 
the strands in the loop does not seem to last for the en- 
tire falling time of the coronal rain. Figure [5] shows that 
after roughly t = 21 minutes it is difficult to observe any 
synchronous oscillation of the strands. This loss of col- 
lective behavior may be apparent, resulting from several 
factors such as, a change in the magnetic field geometry 
of the lines, a change of the morphology of the condensa- 
tions (increasing the difficulty of differentiating between 
the strands), or a change of the plasma conditions in the 
condensation (reheating or further cooling, or a thinning 
of the plasma leading to a change in the opacity in the 
Ca II H line). If the loss of collective behavior is however 
real, it would exclude the option of a kink mode in a sin- 
gle m onolithic loop as it should rema in collective with the 
loop (|Van Doorsselaere et al.ll2QQ8a[ ). For the other two 
options, i.e. kink modes in loop strands and a torsional 
Alfven wave in a single loop, the loss of collective behav- 
ior can be explained in terms of phase-mixing. For the 
former option, variations in the density of the different 
strands leads to a drifting out of phase of the individ- 
ual blob oscillations. For the latter option, if the tor- 
sional modes belong to different magnetic surfaces inside 
the loop, different Alfven speeds are expected, leading in 
turn to a loss of phase. This could explain, in turn, the 
obtained different periods in Table [1] and different phase 
s peeds in Table [p 

iTerradas et al.l (|2QQ8l ) have recently shown analytically 
that a standing kink mode in a loop where siphon flows 
are present will experience a linear shift in phase with po- 
sition along the loop and an asymmetric profile in time of 
the eigenfunctions with respect to the loop's apex. The 
short interval of time during which the loop can be ob- 
served does not allow us to detect any overall phase shift 
for the observed strands. However, due to the relatively 
low speed of the detected flow (< 60 km s~^), we do not 
expect its effect on the phase speed of the waves to be 
significant. For instance, if our oscillations correspond 
to a shifted fundamental mode we would still be able to 
observe the coronal rain oscillating significantly at the 
apex for several periods, which is not observed. 

Figure shows a zoom-in of condensations 1 and 2 
(and 6) while they fall. The x-axis corresponds to the 
projected distance from the apex of the loop (coded in 
color in Figure [5|), and the ^-axis corresponds to the same 
transversal width across the loop as in Figure [5l The 
corresponding times are set in the lower left corner of 



each panel. The 2 minute sequence corresponds approx- 
imately to one period of strands 1 and 2, when these 
exhibit the maximum oscillation amplitude. As can be 
seen in the Figure, both condensations start horizontal 
(time t = 17.25 min), are slanted halfway through the 
oscillation {t = 18.25 min), and end up again horizon- 
tal at the end {t = 19.25 min). Under the assumption 
that the condensations retain their shape during the 2 
minutes of the sequence, so that we are indeed following 
the same plasma parcel, the later means that the con- 
densations do not oscillate up and down uniformly but 
rather exhibit amplitude differences along their lengths. 
As shown in the sketch of Figure [HI in the scenarios of- 
fered by both a kink mode and a torsional Alfven wave 
it is hard to explain such effect only on the basis of, 
respectively, the displacement of the loop axis and the 
twisting of the magnetic field lines. This is due to the 
fact that the observed wavelength of 80 Mm (in case 
of a first harmonic) is much longer than the length of the 
condensations (< 3 Mm), and also to the fact that the 
phase speed of the wave (> 400 km s~^) is much faster 
than the falling speed of the condensation (< 60 km s~^), 
and hence we should expect the whole condensation to 
displace transversally in a uniform way along its length. 
In the case that the length of the condensation is not 
so short when compared to the wavelength (as in the 
sketch of Figure [H we would expect a periodical change 
of the slope of the condensations (with respect to the 
axis of the loop), which is not observed in our case. If 
we have instead a propagating wave, different portions of 
a condensation could oscillate with different amplitudes 
at a given time, but eventually all portions of the con- 
densation should exhibit the same maximum amplitude, 
which is not observed neither. Hence this effect may not 
be caused by the nature of the oscillations. 

We believe the cause for the observed effect in Fig- 
ure [3 to be linked to the cause of the often observed 
separation (and subsequent elongation) of the plasma in 
the condensation. Due to the high density of the con- 
densation it is possible that the plasma beta parame- 
ter is high enough that the plasma moves transversally 
to the axis of the loop, thus allowing also the observed 
separation process. The density and temperature range 
of coronal rain are not well known observationally, but 
since its opacity is large enough to appear bright and 
dark in Ha and Ca II H when observed a bove limb and 
on disk (towards th e limb) respectively (| Antolin et al.l 
12011. in preparation! ), the range of values must be chro- 
mospheric. This is also suppo rted by numerical s i mula- 
tions of catastroph ic cooling (|Muller et al.l I2003L 120041 : 
I Antolin et al.l l2010f ) . In the later the temperature of 
coronal rain is estimated to be as low as 6 x 10^ K 
and its number density to be about 10^^ cm~^. Tak- 
ing a coronal magnetic field inside the loop of 14 G, 
an average of the values calculated in Table [21 we ob- 
tain a rough estimate for the plasma beta parameter of 
/3 = Sirp/B^ SirnkBT/B^ 0.1. Since the loss of 
thermal equilibrium implies high velocity upflows and 
subsequent shocks, it is not unreasonable to consider the 
possibility of the plasma beta parameter becoming high 
enough so that the plasma expands transversally. 3D nu- 
merical simulations of coronal rain formation are needed 
to correctly address this idea, which is the subject of a 
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future work. 

After separation of the initially dense condensation, 
the plasma in coronal rain is observed to elongate into 
strand-like structures, thus probably tracing the internal 
structure of coronal loops. Whether the elongation pro- 
cess is just a result of gravity acting differentially along 
the magnetic topology or if other more sophisticated pro- 
cesses are involved is an interesting question that needs 
to be addressed properly with the help of numerical sim- 
ulations. For instance, since the region below the falling 
blob is expected to be magnetically dominated, accord- 
ing to the previous discussion the layer in between should 
meet the criterium for the onsetting of the magnetic 
Ray leigh- Taylor instability, which would contribute to 
the separation of the blob. This effect is thought to 
be responsibl e for the finder-like structures observed in 
prominences (jBerger e t al. 2010^. 

The causes for the observed oscillations are less clear. 
Reported horizontal oscillations of loo ps are often linked 
directly or indirectly to flares or CMEs (|Nakariakov et al.l 
12009). However, in our observations no flares or ener- 
getic events were reported on that day. A possible cause 
may be interaction with neighboring loops. In Figure [1] 
we have outlined in dotted curves some paths of coronal 
rain marking the presence of other loops. Due to the 
projection effect it is hard to say whether these loops are 
indeed close-by and whether any interaction really oc- 
curs. It is interesting to note however that the coronal 
rain in these loops occurs prior to the coronal rain in 
our loop, and that the paths seem to intersect roughly 
halfway through the visible leg of our loop. If this is not 
a projection effect, it is possible that the coronal rain 
perturbs our loop thus producing the oscillations. The 
perturbation would have a maximum amplitude at the 
crossing region, thus explaining why we observe a max- 
imum amplitude for the oscillations halfway along the 
leg. 

Another possibility can be a scenario in which the in- 
ertia of the condensations conforming coronal rain is not 
negligible, thus affecting the stability of the entire loop 
(and hence triggering themselves oscillations in the loop). 
Making a rough analogy to a hose with water gushing in, 
it is natural to expect that there might be a limit for the 
quantity and for the velocity of the plasma above which 
the stability of the magnetic field structure is compro- 
mised and the loop oscillates. This scenario is worth 
investigating through numerical simulations, and will be 
the subject of a subsequent paper. 

An additional possible scenario for triggering transver- 
sal oscillations in loops is one in which the waves 
are linked to an unobserved energetic event such 
as magnetic reconnect ion in the lower atmosphere, 
triggered by convection. 3D Numerical simulations 
have shown that magnetic reconnection processes 
in the chromosphere (in a solar atmosphere pow- 
ered by convection) can easily generate energetic 
events such as spicules and o t her je t- like phenomena 
iMartmez-Svkora^eiaU I2009L 120101 : iHeggland et~all 
[2011, in prepara tioiJ ) . It is thus reasonable to con- 
sider reconnection as a possible candidate to generate 
transverse perturbations of flux tubes. On the other 
hand, swirl events have been observed in the photos phere 
([Wedemever-Bohm fc Rouppe van der Voortl[2009h . Nu- 
merical simulations have shown that torsional motions at 



photospheric level can gen erate different ki nds of mode s 
in flux tubes iFedun et al.l ( Submitted on 2 7 Jan 2QTT| ): 
iShelvag et al.l ([Submitted on 31 Jan 20111) and in- 
put enough amounts of energy for coronal h eating 
(jKudoh fc ShibatalfToollAntolin fc Shib ata 20101. 

Since we expect a certain degree of twisting and braid- 
ing for the strands in the observed loop, a last possible 
scenario we can think of is one in which the observed 
oscillatory events are not due to waves, but result from 
the internal complex topology of the loop. The conden- 
sations would then fall down the braided 'helical' strands 
which would make it seem they are oscillating. We con- 
sider however this scenario to be unlikely. Since we ob- 
serve the oscillations over several strands, and mostly 
in-phase, the twist in the loop would need to exist over 
most of the loop's width. Also, we would need the tube 
to be twisted along its entire length in order for it to 
be in equilibrium in the corona. However, the twist we 
observe seems to have 2 nodes, one at the apex and one 
towards the footpoint, and its amplitude increases in be- 
tween. Furthermore, the stability of such a configuration 
has to be maintained over the time interval in which the 
oscillatory events are observed, which is longer than 20 
minutes. 

4. CONCLUSIONS 

We have analyzed transversal oscillations of a loop that 
are put in evidence by coronal rain falling in the loop. 
The coronal rain is observed to fall down from the apex, 
roughly 25 ± 5 Mm above the surface. The condensa- 
tions composing coronal rain are observed to separate 
and elongate as they fall down, exhibiting large distribu- 
tions of velocities but rather concentrated accelerations 
around a much lower value than that of the average ef- 
fective solar gravity along the loop, thus implying the 
presence of a different force, probably of magnetic ori- 
gin. The obtained deceleration can be the result of wave 
pressure from transverse magnetohydrodynamic waves at 
coronal heights. 

As the condensations fall, they elongate into strand- 
like structures that are observed to oscillate in-phase 
transversally to the axis of the loop with periods that 
are similar to those normally observed in prominences. 
The amplitudes of the oscillations are observed to vary 
significantly with respect to the position along the loop, 
having a maximum at roughly halfway through one leg 
and minimums at both the apex and towards the foot- 
point of the leg. We have interpreted this result as a 
signature of a first harmonic of a standing transverse 
magnetohydrodynamic wave in the loop, although an up- 
ward propagating wave is also a possible, but less likely, 
scenario. This interpretation implies a wavelength equal 
to the loop's length, 80 ± 15 Mm. 

Since this active region loop exhibits a catastrophic 
cooling event we expect the internal density to be sig- 
nificantly higher than that of the external corona for an 
interval of time long enough to create a waveguide along 
the loop. The obtained phase speeds of the waves are 
between 400 km s~^ and 1000 km s~^, implying either a 
fast (horizontal) kink mode or a torsional Alfven mode. 
The wide distribution of speeds may be due to the pos- 
sible uncertainties in the measurements, given the short 
time in which the loop can be observed. On the other 
hand, if the distribution in the phase speeds is real, it im- 
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plies a loss of collective behavior which can be explained 
in terms of phase-mixing, present both in a scenario in 
which each strand has its own kink mode, a nd in the sce- 
nario of a torsional Alfven wave. Recently, iHarris et alJ 
()2Qllf ) reported collective behavior in transverse oscil- 
lations of an ensemble of prominence threads. There 
also the collective nature is lost with time as the threads 
oscillate with slightly different periods and amplitudes. 
The average coronal magnetic field inside the loop is esti- 
mated to be between 8 ± 2 G and 22 ± 7 G, in agreement 
with other estimates of coronal magnetic fields in ac- 
tive region loops through co ronal seismology techniques 
(|Nakariakov fc Ofmanll2QQTl ). 

The strand-like condensations do not exhibit a uniform 
oscillation amplitude along their lengths, which we be- 
lieve is not caused by the oscillation, but by the physical 
conditions inside the loop. A rough estimate indicates 
an average plasma- /3 parameter of 0.1 in the condensa- 
tion due to the high densities and strong shocks that 
are normally created by the catastrophic cooling mecha- 
nism. This allows us to interpret the result as evidence 
of strong gas pressure forces relative to magnetic forces 
inside the loop, which would explain the often observed 
initial separation of the condensations and which can also 
be linked to the deceleration of the plasma. 

Coronal rain has been previously shown to be deeply 
linked with the coronal heating mechanism in the loop. 
Here we have shown the potential it can play in under- 
standing the magnetic field topology of the solar corona 
by tracing the internal substructures of loops, marking 
the internal forces at play, exposing wave-like phenom- 
ena and thus allow the measurement of the coronal mag- 
netic field strength by means of coronal seismology. We 
have pointed to several important problems that need 
to be addressed in future work. Namely, by means of 
3D numerical simulations, investigation of the possible 
processes allowing the separation and elongation of the 
condensations in the coronal rain, investigate the inter- 
nal physical conditions of such condensations created 
through the catastrophic cooling mechanism by means 
of a proper radiative transfer model of the atmosphere, 
and investigate coronal rain being itself a possible cause 
of transverse magnetohydrodynamic oscillations of loops. 
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